Forward-Backward Multiplicity Correlations in y/s^N = 200 GeV Au+Au Collisions 
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Forward-backward correlations of charged-particle multiplicities in symmetric bins in pseudo- 
rapidity are studied in order to gain insight into the underlying correlation structure of particle 
production in Au+Au collisions. The PHOBOS detector is used to measure integrated multiplici- 
ties in bins centered at n, defined within \rj\ < 3, and covering intervals A-q. The variance oq of a 
suitably defined forward-backward asymmetry variable C is calculated as a function of r/, An, and 
centrality. It is found to be sensitive to short range correlations, and the concept of "clustering" is 
used to interpret comparisons to phenomenological models. 
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At the top energies reached for Au+Au collisions at 
the Relativistic Heavy Ion Collider (^/snn — 200 GeV), 
most of the energy is carried by particles with large 
longitudinal momenta. In general, these momenta are 
best expressed with the rapidity variable (y = 4 In gjr^ ) 
or its near equivalent, the pseudorapidity variable (rj — 
— ln(tan(0/2))). In strongly interacting systems, it is of- 
ten thought that correlations between two particles are 
mainly "short-range" in rapidity. Indeed, short range 
correlations have been observed in pv and pp collisions 
over a wide range of beam energies |l| via two-particle 
correlation measurements. However, they have never 
been shown to be the only source of correlation in mul- 
tiparticle production. 

Indeed, single-particle rapidity (or pseudorapidity) dis- 
tributions for inclusive charged particles produced in 
heavy ion collisions reveal a distinct "trapezoidal" struc- 
ture stretching across the full rapidity range available 
in these reactions. More importantly, there is no ev- 
idence for an extended boost invariance, which would 
imply independent emission from different rapidity re- 
gions. Rather, two non-trivial effects are visible in the 
centrality dependence of particle production, which are 
apparently long-range in rapidity. The first is the fact 
that integrating over the full phase space reveals that 
particle production depends linearly on the number of 



participants |2j. The second is that this occurs despite 
a significant change in the shape of the pseudorapidity 
dependence as a function of centrality, which happens to 
be collision-energy independent in the forward region [j| . 

All of this suggests that charged-particle production 
is highly correlated over large regions of rapidity, which 
naturally begs the question of the underlying structure 
of the single-particle distributions. In this paper, we dis- 
cuss how to take first steps in this direction via the study 
of forward-backward multiplicity correlations. By this is 
meant the event-by-event comparison of the integrated 
multiplicity Np in a bin defined in the forward (n > 0) 
region, centered at rj with pseudorapidity interval Arj, 
with the multiplicity Nb measured in an identical bin de- 
fined in the backward hemisphere, centered at —77. With 
these definitions, one can construct the event-wise ob- 
servable C — (Np — Nb)/^Nf + Nb, and measure the 
variance for a set of events with nominally similar 
characteristics (e.g. collision centrality). 

The C variable is chosen to have particular sensitivity 
to various types of long and short range correlations. An 
"intrinsic" long-range correlation in the emission of 
particles into the forward and backward hemisphere from 
a single source would give Np—Ns = with a substantial 
value of Np + Nb, forcing C to 0. However, if the parti- 
cle sources tended to produce into the forward or back- 
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ward region, such that the partitioning was binomial, this 
would lead to = 1, since cr 2 (N F ~ N B ) = N F + N B 
in that case. Direct studies of (Np) as a function of 
Nb were used by UA5 for pp collisions to conclude that 
there are no large intrinsic correlations in particle pro- 
duction. And yet, a non-trivial long range correlation 
between the hemispheres persists, even when excluding 
the central two units of pseudorapidity Q, 0] . 

Short-range correlations would arise if objects emit- 
ted into either hemisphere break into (k) particles on 
average, with a variance of o~\, each of which remains 
close in rapidity (e.g. due to isotropic emission). Such 
intrinsic short-range correlations have in fact been mea- 
sured in pp and pp experiments by direct construction 
of the two-particle correlation function in 77. In a series 
of papers, the UA5 collaboration explored the hypoth- 
esis that particle emission was dominated by the decay 
of "clusters" locally in rapidity space Q, . By varying 
A?7 and measuring the ratio Aap / (Np) as a function of 
Np + Nb , they found an "effective" cluster multiplicity 
(k e ff = (k) + a\f (k)) of approximately 2 charged parti- 
cles. This exceeds the multiplicity expected from a res- 
onance gas using masses up to 1.5 GeV, which gives an 
average multiplicity per particle of about (k) ~ 1.5 
suggesting that clusters may have a variety of dynamical 
sources (although, to be sure of this, resonances over a 
wider mass range than in the previous calculations should 
be taken into account). 

The decay of clusters has a particular effect on Oq. 
Consider the idealized case, where a cluster decays into 
exactly k particles, and all of the N particles are ran- 
domly distributed into identical fy-bins in the forward or 
backward directions. In this case, the underlying fluctu- 
ations are those associated with N/k objects rather than 
N independent particles, and C — > \f~k~C. In other words, 
is linear with the cluster multiplicity for this simple 
case. In a more realistic physical situation, particles from 
each cluster may fall outside the chosen 77-bin, or even in 
the opposite hemisphere. This results in a non-trivial 
modification of the measured value of k, but one which 
is still related directly to k, or rather k e ff (which incor- 
porates the effect of the multiplicity per cluster having its 
own distribution) . This relationship can be explored with 
the use of simulations and thus used to extract the effec- 
tive cluster multiplicity with measured forward-backward 
correlations. These studies show that forward-backward 
correlations have the surprising ability to measure a fun- 
damental local property of particle production. 

The data analyzed here were taken with the PHOBOS 
detector @ during Runs 2 and 4, in 2001 and 2004, re- 
spectively. The pseudorapidity acceptance was restricted 
to that of the "Octagon" detector, which is a tube of sili- 
con sensors covering \rj\ < 3 and full azimuth except for a 
region near midrapidity. To simplify the analysis for dif- 
ferent values of rj and Arj, only the regions of the detec- 
tor with complete rapidity coverage were kept, restrict- 
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FIG. 1: Number of hits vs. total deposited energy (in keV) 
measured —0.5 < 77 < (in the Octagon detector) for 200 GeV 
Au+Au collisions. The solid line is a function, described in 
the text, fit to the mean number of hits in each energy bin. 



ing the total azimuthal acceptance to A<f> — it, in four 
45-degree wedges • The multiplicity in each bin is esti- 
mated event-by-event by summing up the angle-corrected 
deposited energy of all hits and then dividing by the aver- 
age energy per particle ||. The average is calculated by 
a fit to the distribution of (Nhits) vs. the total deposited 
energy (E), an example of which is shown in Fig.^ with 
the function (N hit ) = N max (l - e - E l E ^). The aver- 
age energy per particle is given by E max /N max . This 
removes nearly all occupancy-related effects. Beyond the 
usual lower threshold to define a hit in the silicon, an in- 
dependent upper bound of the deposited energy per hit is 
applied to reduce the effect of low- momentum secondaries 
on the fluctuations of the single-particle dE/dx distribu- 
tions, while keeping 97% of the primaries. The average 
(C), calculated as a function of 77, centrality, and event 
vertex, is subtracted event-by-event to correct for gaps 
in the Octagon. This correction was studied in detail 
with simulations, and was found to leave the fluctuations 
unaffected. 

To provide information that can be compared directly 
to models, a procedure was developed to estimate and 
remove the detector effects from the raw measured value 
of a 2 c {o% raw ) by using Monte Carlo (MC) simulations 
of the PHOBOS apparatus. The basic idea is to assume 
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where a^ et is the contribu- 
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tion from detector effects, and use the MC to subtract 
it on average, leaving behind only the correlations be- 
tween primary particles. It is found that there are sev- 
eral sources which contribute differently as a function of 
77 and combine in quadrature to a nearly-constant value 
over the pseudorapidity range covered by the Octagon |]| , 
as illustrated in Fig. [21 Gaps between the silicon sen- 
sors tend to be most important near mid-rapidity and 
decrease rapidly with angle as the particle density de- 
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FIG. 2: Contributions to oq from sources other than primary 
charged particles. 




FIG. 3: Example of the analysis method, showing the raw 
values of erf. (erf; rau ,), the contribution from detector effects 
(a| et , as estimated using modified HIJING), and the sub- 
tracted result (erf.) for bins of A77 = 0.5 as a function of r\ in 
the Run 4 data set. Note that these values are calculated for 
the half-azimuth (Ac/> — n) acceptance, which is corrected-for 
in the final result. 

creases. Similarly, larger-angle particles tend to be slower 
than those at forward angles, leading to larger fluctua- 
tions in the energy deposition in the single layer of sili- 
con. At the same time, the number of secondary parti- 
cles emitted from primaries interacting in the beryllium 
beampipe rises quickly with increasing rapidity. This 
leads to increased fluctuations, even though some of the 
background is removed with a properly-chosen lower cut 
on the angle-corrected dE/dx of each hit. We have also 
studied fluctuations from the momentum (or j3) distri- 
butions of the particles incident on the detector. In this 
case, the contribution to erf. decreases with increasing 
rapidity. 

The detector effects are corrected on average by calcu- 



FIG. 4: Calibration of reconstruction method with MC sim- 
ulations, as a function of r\ with A77 = 0.5. Solid points rep- 
resent MC truth information on erf., while the open symbols 
represent the outcome of the reconstruction procedure. 

rating erf. using a modified HIJING simulation with all 
intrinsic correlations destroyed by flipping the sign of r\ 
at random, effectively breaking any rapidity-dependent 
multi-particle correlations and forcing crfi = 1. This 
allows a direct estimation of o\ et . An example of this 
is shown in Fig. which shows the raw values of erf, 
( a Craw)' the contribution from detector effects (o^ et , as 
estimated from modified HIJING), and the subtracted 
result (erf.) for bins of A77 = 0.5 as a function of r\ in the 
Run 4 data set. There is a small residual correlation of 
o-det with °c> parametrized as a\ et = o\ et raw {\-a{a 2 c - 
1)), which is removed by estimation of a as a function of 
77, A77 and centrality. Finally, we correct for using only 
half-azimuth acceptance by erf. — > 2 (erf, — 1) + 1, a for- 
mula obtained using MC simulations, assuming that the 
limited acceptance only reduces the probability of mea- 
suring all of the particles in a typical cluster. 

Systematic uncertainties were calculated by varying 
several variables involved the estimation of erf, (data 
sets, event generators, and dE/dx cuts) and found to 
be around Acrf, ~ 0.1. The bin-to-bin variation of the 
systematic error calculation was adjusted to reduce fluc- 
tuations from the error determination procedure. 

After correcting for detector effects, the results on erf, 
can be directly compared with model calculations based 
on charged primary particles. We have focused mainly 
on HIJING @ and AMPT [HI (which includes partonic 
and hadronic transport codes), both of which have been 
used to describe various features of heavy ion collisions 
at RHIC energies. Fig. 01 shows the results of the cor- 
rection procedure described above on three simulations. 
Two of them are based on standard versions of HIJING 
and AMPT, and a third on the modified HIJING, as de- 
scribed previously. One sees that in all three cases, the 
reconstructed values of Oq, using the same tuning in all 
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cases, match the Oq extracted from the primary charged 
particles within statistical precision. 

The first set of results concerns the r\ dependence of 
Uq, for forward and backward bins that are A77 = 0.5 
units wide, as shown in the upper panel of Fig.0 In this 
case, HIJING and AMPT already show some differences 
related to their underlying physics scenarios. For periph- 
eral (40-60%) events, both models have a similar magni- 
tude and a monotonically-rising rj dependence. Central 
events show a substantial difference extending over most 
of the rapidity range, with AMPT showing a systemati- 
cally smaller value of Oq . This may be due to the initially 
produced clusters being smeared out in rapidity space by 
the hadronic rescattering stage. Also, it is observed that 
both data and MC show a\, ~ 1 at r\ = 0. This is 
suggestive of the effect predicted by Jeon and Shi [llj . 
that the formation of a QGP near mid-rapidity should 
break up any sort of cluster structure seen in pp and thus 
lead to a reduction in er^. However, it can also be ex- 
plained by the fact that clusters produced at rj = tend 
to emit particles into both the Np and Nb side, inducing 
an "intrinsic" long-range correlation that decreases <Jq. 
Simulations show that this effect is not modified by the 
precise values of (k) and erf., for constant fc e //. In any 
case, the data are systematically higher than the model 
calculations for both peripheral and central events. 

The next set of results is the A77 dependence of cTq, 
with a fixed bin position 77 = 2.0, as shown in the lower 
panel of Fig. [SJ One sees that for both peripheral 
and central data, Oq rises monotonically with increas- 
ing pseudorapidity interval. This can be explained as an 
acceptance effect in the context of cluster emission: by 
increasing A77, one increases the probability of observing 
more than one particle emitted from a single cluster in 
either Np or Nb- Clearly, the rate of change with A77 
should reflect the full cluster distribution (both (k) and 
crfe), but these studies are not yet sufficiently precise to 
determine the detailed properties of clusters produced in 
RHIC collisions. It is striking that the peripheral data 
has already reached o-q ~ 3 for the largest A77 while this 
number is closer to 2 for central data. Finally, it is in- 
teresting that neither HIJING nor AMPT can explain 
both the centrality and Ar] dependence simultaneously. 
AMPT never agrees with the data in magnitude, but at 
least predicts a larger <Jq in peripheral events than in 
central events. This should be contrasted with HIJING, 
which reproduces the central data but has no centrality 
dependence at all, at odds with the experimental data. 

In conclusion, measurements of forward-backward fluc- 
tuations of charged-particles produced in Au+Au colli- 



sions at y?SNN = 200 GeV may provide insight into the 
structure of long and short range correlations in pseudo- 
rapidity space. Our data for 200 GeV Au+Au collisions 
are now fully corrected for detector and background ef- 
fects, so direct comparisons can be made to phenomeno- 
logical models. We see significant short-range correla- 
tions at all centralities and pseudorapidities, instead of 
just at mid-rapidity. There is a non-trivial centrality and 
rapidity dependence of these correlations, in both 77 and 
A77. Finally, neither HIJING nor AMPT reproduces all 
of the main qualitative features, but the way in which 
they fail to do so may well provide information on the 
underlying physics. In particular, more theoretical at- 
tention should be paid to the properties of "clusters" re- 
quired to explain our data. As mentioned above, Jeon 
and Shi have proposed that QGP formation would mod- 
ify the measured properties of clusters ^lj. The data 
shown here should provide means to study such effects, 
or set upper limits on their occurrence. 
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FIG. 5: (upper) cr£ for a fixed Ar? = 0.5 as a function of 77. (lower) 0% as a function of A77 for a fixed bin center at 77 = 2.0. 
Statistical errors are shown as vertical bars, while systematic errors are shown as brackets. 



